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Abstract 

We review lattice QCD results for glueballs (including a discussion of 
mixing with scalar mesons), hybrid mesons and other exotic states (such 
as BsBs molecules). 

1 Introduction 

Quantum Chromodynamics has emerged as thc unique theory to describe hadronic 
physics. It is formulated in terms of gluonic and quark fields. The only free 
parameters are the scale of the couphng (usually called Aqcd) and the quark 
masses defined at some conventional energy scale. 

Where large momentům transfers occur, the effective coupling becomes weak 
and a perturbative treatment is valid: in this domain the theory has been tested 
directly by experiment. However, because the effective coupling is weak for these 
processes that can be described by perturbation theory they are necessarily not 
the dominant hadronic processes. A typical hadronic process will involve small 
momentům transfers and so has to be treated non-perturbatively. 

In this non-pcrturbative régime, the description of hadrons is quite far re- 
moved from the description of the gluonic and quark fields in the QCD La- 
grangian. Because only colour-singlet states survive, the hadrons are all com- 
posites of quarks and gluons. One example emphasises this: the nucleon has 
a mass which is very much greater than the sum of the quark masses of the 
three valence quarks comprising it. This extra mass comes from the gluonic 
interactions of QCD. Another way to view this is that the naiVe quark model 
is a useful phenomenological tool but has constituent quarks with masses much 
greater than the QCD masses (ie masses as defined in the Lagrangian). It is 
important to understand why this is approximately what QCD requires and to 
find where QCD departs from the naivě quark model. 

One way to characterise the manner in which QCD goes beyond thc naivě 
quark model is through the concept of exotic states. Here exotic is taken to mean 
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'not includcd in the naivě quark model'. In ordcr to discuss exotic states, we 
necd to summarise what the naiVe quark model contains. Basically the degrees 
of frcedom are the valence quarks (ie quark-antiquark for a meson and 3 quarks 
for a baryon) with masses and intcractions given by some effective interaction. 
Tho conscqucnccs of this are that only ccrtain J^*^ values will exist and that 
the munber of states with different quark flavours is specified. So, concentrating 
on mesons madc of the three flavours of light quarks (u, d, s), onc expccts a 
nonet of mesons with the flavours (ud, du, uu±dd, šs, us, ds, šu, šd). This is 
indeed what is found for vector mesons (p, w, </>, K*). It is also possible within 
the quark model for the flavour-singlet states {uu + dd, šs) to mix, as found for 
the pseudoscalar mesons. What would be exotic is for a tenth statě to exist. 
For mesons with orbital angular momentům L between the quark and antiquark 
the allowed J^^ values are shown below. Thus spin-parity combinations such 
as O , O"*"", 2^ are termed spin-exotic since they cannot be made from 
a quark plus antiquark alone. 

It has been a considerable challenge to build a machinery that allows non- 
perturbative calculations in QCD with all systcmatic errors determined. The 
most controUed approach to non-perturbative QCD is via lattice techniques in 
which space-time is discretized and time is taken as Euclidean. The functional 
integrál is then cvaluatcd numerically using Monte Carlo techniques. 

Lattice QCD nccds as input the quark masses and an overall scale (conven- 
tionally given by Aqcd)- Then any Green function can be evaluated by taking 
an average of suitable combinations of the lattice fields in the vacuum samples. 
This allows masses to be studied easily and matrix elements (particularly those 
of weak or electromagnetic currents) can be extracted straightforwardly. Un- 
like experiment, lattice QCD can vary the quark masses and can also explore 
different boundary conditions and sources. This allows a wide range of studies 
which can be used to diagnose the health of phenomenological models as well 
as casting light on cxperimental data. 

One limitation of the lattice approach to QCD is in exploring hadronic decays 
because the lattice, using Euclidean time, has no concept of asymptotic states. 
One feasible stratégy is to evaluate the mixing between states of the same energy 
- SO giving some Information on on-shell hadronic decay amplitudes. 

There is an interesting theoretical world in which the quark degrees of free- 
dom are removed from QCD, leaving pure gluo-dynamics. This is also known 
as pure Yang-Mills theory. It is a self-consistent theory which has the fuU non- 
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perturbative gluonic interaction. It turns out that this gluonic interaction does 
produce the salient features of QCD: asymptotic freedom, confinement, etc. It 
is of interest to explore the spectrum in this casc: the states are called ghieballs. 

It is also of interest to consider the propagation of quarks in this gluonic 
theory. The quarks are treated as in the Dirac equation and they propagate 
through the ghionic ground statc. This approach is known as the qucnchcd 
approximation. Again this turns out to be a very useful approximation: chiral 
symmetry breaking occurs for example. This quenched approximation is in 
contrast to the fuU quantum field theory of QCD where there would be quark 
loop effects in the ground statě also. Thus in the quenched approximation there 
will be inconsistencies: the theory is not unitary. However, for heavy quarks 
it will be a good approximation since heavy quark loops are suppressed and it 
may be adequate to describe some features of lighter quarks. Moreover, many 
phenomenological models are appropriate to the quenched case and so can be 
compared with quenched QCD. 

2 Glueballs and scalar mesons 
2.1 Glueballs in quenched QCD 

Glueballs are dcfincd to be hadronic states madc primarily from gluons. The 
full non-perturbative gluonic interaction is included in quenched QCD. A study 
of the glueball spectrum in quenched QCD is thus of great value. This will allow 
experimental scarchcs to be guidcd as well as providing calibration for models 
of glueballs. A non-zero glueball mass in quenched QCD is the "mass-gap" of 
QCD. To prove this rigourously is one of the major challenges of our times. 
Herc WC will explore tlic situation using computational tcchniqucs. 

In lattice studies, dimensionless ratios of quantities are obtained. To explore 
the glueball masses m, it is appropriate to combine them with another very 
accurately measiircd quantity to have a dimensionless observable. Since the 
potential between static quarks is very accurately measured from the lattice, 
it is now conventional [1] to use ro for this comparison. Here ro is implicitly 
defined by r^dV{r)/dr = 1.65 at r ro where V{r) is the potential energy 
between static quarks which is easy to determine accurately on the lattice. 
Conventionally ro ~ 0.5 fm. 

Thcorctical analysis indicatcs that for Wilson's discretisation of the gauge 
fields in the quenched approximation, the dimensionless ratio mro will diífer 
from the continuum limit value by corrections of order a?. Thus in fig. 1 the mass 
of the J^^=0++ glueball is plottcd versus the lattice spacing a . The straight 
line then shows the continuum limit obtained by extrapolating to a = 0. As can 
be seen, there is essentially no need for data at even smaller o-values to further 
fix the continuum value. The value shown corrcsponds to m(0'''"'')ro = 4.33(5). 
Since several lattice groups [2, 3, 4, 5] have measured these quantities, it is 
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Figuře 1: The value of niass of the J = 0++ glueball statě from quenched data 
[Np = 0)[2, 3, 4, 5] in units of ro where vq w 0.5 fm. The straiglit line shows 
a fit describing the approach to the continuum limit as a ^ 0. Results [6, 7, 9] 
for the lightest scalar meson with Np = 2 flavours of sea quarks are also shown. 



reassuring to see that the purely lattice observables are in excellent agreement. 
The publicised difference of quoted m(0++) from UKQCD [4] and GFll [5] 
comes cntirely from relating quenched lattice measurements to vahies in GeV. 

In the quenched approximation, diíferent hadronic observables differ from 
experiment by factors of up to 10%. Thus using one quantity or another to 
set the scale, givcs an overall systematic error. Herc the scale is set by taking 
the conventional value of the string tension (determined from potential models 
and from hadronic lattice studies), ^/a = 0.44 GeV, which then corresponds to 
= 373 MeV (or ro = 0.53 fm). An overall systematic error of 10% is then 
to be included to any extracted mass. This yields to(0"'""'") = 1611(30)(160) 
MeV where the second error is the systematic scale error. Note that this is 
the glueball mass in the quenched approximation - in the real world significant 
mixing with qq states etc may modify this value substantially, as we discuss 
below. 

In the Wilson approach, the next lightest glueballs are [3, 4] the tensor 
m(2++)ro = 6.0(6) (resulting in m(2++) = 2232(220)(220) MeV) and the pseu- 
doscalar m(0 '")ro = 6.0(1.0). Although the Wilson discretisation provides 
a definitivě study of the lightest (0++) glueball in the continuum limit, other 
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Figuře 2: The continuum glueball spectrum[10]. 



methods are competitive for the determination of the mass of heavier glueballs. 
Námely, using an improved gauge discretisation which has even smaher discreti- 
sation errors than the dependence of the Wilson discretisation, so aUowing a 
relatively coarse lattice spacing a to be used. To extract mass values, one has to 
explore the time dependence of correlators and for this reason, it is optimum to 
use a relatively small time lattice spacing. Thus an asymmetric lattice spacing is 
most appropriate. The results [10] are shown in fig. 2 and for low lying states are 
that m(0++)ro = 4.21(11)(4), m(2++)ro = 5.85(2)(6), m{0-+)ro = 6.33(7)(6) 
and m{l^ )ro = 7.18(4)(7). It will be very diflicult to identity experimentally 
states corresponding to these heavier glueballs since the spectrum is rich in qq 
states of those quantum numbers at those mass values and there will thus be 
considerable mixing. 

One signál of great interest would be a glueball with J^'~^ not allowed for 
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qq - a spin-exotic glueball or oddball - since it would not mix with qq states. 
These states are found [3, 4, 10] to be high lying: considerably above 2m(0++). 
Thus thcy are likely to be in a region around 4 GeV where it is very difRcult to 
separate states unambiguously by experiment. 

As well as the mass of a glueball, it is possible to study their physical size. 
In principle this is detcrmined by measuring the matrix element {G\J\G) where 
J is some local current which couples to the glueball. Since glueballs have 
no flavour, the energy-momentum tensor is the most appropriate choice. A 
preliminary study has been made, albeit with large systcmatic crrors [12] and 
finds a rádius 0.9 ± 0.3 fm. This approach should be contrasted with what 
has come to be called the Bethe-Saltpeter wavefunction which is obtained from 
(G|-L|0) where L is the lattice operátor used to crcate a glueball statě from the 
vacuum. Within a lattice calculation, it is easy to measure the dependence of 
this overlap on the spatial extent of L [11], but difíicult to interpret the result. 

Related Information can be obtained from a study of the gluelump: the 
statě with one static colour source in the octet representation with a gluonic 
field making it a colour singlet. This would be of physical relevance should a 
massive gluino exist: it would be the glueballino, a gluino-gluon bound statě. 
The spectrum [13] and spatial distribution [14] have been studied. 

Another topič which is of mainly theoretical interest is the glueball mass 
(as a dimensionless ratio to the string tension) as the number of colours Nc is 
varied from 3. The SU(2) Yang-Mills theory has often been studied, especially 
as it is computationally simpler. Recently a study of SU(N) for N=4 and 5 has 
been made. The summary [22] is that N=oo is relatively close to N=3. This 
has theoretical implications since the N=oo theory is formally simpler. For a 
comparison of these lattice results with ADS supergravity see ref. [23] . 

2.2 Scalar mesons in quenched QCD 

In quenched QCD the flavour singlet (/o) and non-singlet (ao) scalar mesons are 
degeneráte. In fuU QCD this degeneracy is split by disconnected quark diagrams 
but these are omitted from the quenched approximation. This same fcaturc of 
the quenched approximation implies that the 7] meson is wrongly treated - it 
will be degeneráte with the tt. This implies that the scalar meson propagation 
can have the wrong sign [18] because the rjn intermediate statě is mistreated 
(once quark loops are allowed in the vacuum then this anomaly is removed). 
For light quarks of mass corresponding to the strange quark or heavier, it is 
expected that this anomaly is relatively unimportant. Thus the measurement 
of the mass of the qq scalar meson can be particularly unreliable in the quenched 
approximation. 

Even though the mixing of the glueball and qq states is not implemented 
in the quenched approximation, one can determine the mixing matrix element. 
This can then be used to estimate the result of the mixing by hand (by using 
a mass matrix for example). On a rather coarse lattice (a~^ w 1.2 GeV), two 
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groups have attempted to measure this inixing [16, 7]. Their results expressed 
as the mixing for two degeneráte quarks of mass around the strange quark mass 
are similar, námely E « 0.3 GcV [16] and 0.5 GeV [7]. This is a relatively large 
mixing (if the glueball and scalar meson states were degeneráte they would be 
split by ±E). 

An exploratory attempt to extrapolate this mixing to the continuum [16] 
gave a very small mixing of 61(45) MeV, while the other determination [7] uses 
clover improvement so order a effects in the extrapolation to the continuum 
are supprcssed and onc would not cxpcct a significant dccrcasc in going to the 
continuum limit. What this discussion shows is that precision studies of the 
mixing on a quenched lattice have not yet been achieved. Furthermore the 
probicms with the scalar meson propagation in the quenched approximation 
discussed above also limit progress. 

As well as this mixing of the glueball with qq states, there will be mixing with 
qqqq states which will be responsible for the hadronic decays. A hrst attempt 
to study this [17] at a coarse lattice spacing yields an estimated width for decay 
to two pseudoscalar mesons from the scalar glueball of order 100 MeV. A more 
realistic study would involve taking account of mixing with the nň and sš scalar 
mesons as well. 

2.3 Scalar mesons in full QCD 

It is now feasible to explore the flavour-singlet scalar meson spectrum including 
the quark loops in the vacuum, ie in full QCD. From dynamical fermion studies 
with Nf = 2, one can determine the flavour singlct and non-singlet mass spec- 
trum. What is found [7, 9] is that the lightcst flavour-singlet scalar meson (/o) 
is lighter than the lightest flavour non-singlet (ao). 

The interpretation of this study is hampered by the same issue that hampers 
the interpretation of cxpcrimcntal data, námely, the mass eigenstates are not 
distinguished as 'glueball' or as 'quark-antiquark'. What one can do is explore 
the output spectrum and deduce what mixing might have occurred. To give an 
example, where we restrict here to Nf = 2 flavours of degeneráte quarks, the 
/o masses will be mo and toq' where the latter is the first excited statě and the 
flavour non-singlet oq mass will be mi. Results for mo are given in fig. 1. Then 
one would expect in a simple 2x2 mixing scenario (ie glueball and qq meson) a 
mass matrix 

f mo E \ 
\ E mi J 

where uig is the glueball mass and E the mixing matrix element. This will have 
two mass eigenstates which can be identified with mo and mo' so determining 

the two frče parametcrs in the matrix. This approach explains what is going on 
but obtains two numbers with two parameters, so there is no cross-check. 
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One can directly address the issue of the mass of the lightest scalar singlet 
meson from the lattice with Nf = 2. It is advantageous to use as full a basis 
of lattice opcrators as possiblc, including Wilson loops and quark-antiquark 
loops. Including the latter can in practice lead to a lower value the ground 
statě scalar meson mass - see refs.[8, 9]. Most studies have shown no significant 
change of the scalar glueball mass as dynamical quarks are includcd [6] . Howcvcr 
the larger lattice spacing result [7] shows a significant reduction in the lightest 
scalar mass, as shown in fig. 1. Before concluding that this implies a lower 
scalar mass in the continuum limit, one needs to chcck whcthcr an cnhanced 
order correction might be present. The origin of the large coeíficient of 
in quenched glueball studies is usually ascribed to the presence of a critical 
point in the fundamcntal-adjoint coupling plane which is closc by in the usual 
Wilson approach with zero adjoint coupling. The extent to which this will be 
enhanced/reduced when dynamical quarks are introduced is not clear. Studies 
using the same approach at a finer lattice spacing [8, 9] do suggest that this 
large order effect is significant for dynamical quarks, but studies even nearer 
to the continuum or with improved actions are needed to resolve this fully. 

A further complication is that as the quark mass is reduced towards the phys- 
ical light quark mass, the decay to nn becomes energetically allowed. The study 
of unstable particles is a difRcult problém in a Euclidean time formalism [19]. 
We return to this topič later. 

2.4 Experimental evidence for scalar mesons 

In full QCD, for the favour-singlct states of any givcn J^*^, there will be mixing 
between the sš statě, the uu + dd statě and the giuoball as well as with multi- 
meson channels. It may indeed turn out that no scalar meson in the physical 
spectrum is primarily a glueball - all states are mixtures of glue, qq, qqqq, etc. 

To help with understanding the experimental situation [20] , we first discuss 
the flavour non-singlet states, the oq with isospin 1. The observed states are 
at 980 and 1450 MeV. The lighter statě has dynamics which appears to be 
dosely associated with the KK thrcshold. The hcavicr statě is not yet very 
well established but seems to be a candidate for a statě mostly comprised of qq, 
while the lighter statě would be qqqq. 

The flavour singlet states (/o) are more numerous. There is a very broad 
enhancement in the tttt S-wave phase shift around 700 MeV (sometimes called 
the a), there is a statě near the KK threshold at 980 MeV and there are more 
states at 1370, 1500 and 1710 MeV. Again assuming that the statě at 980 MeV 
is predominantly qqqq, this suggests that the three states in the 1300-1750 MeV 
energy range are admixtures of the glueball, uu + dď and sš. The fact that there 
are indeed three states in this energy region close to the quenched glueball mass 
of 1600 MeV is the strongest evidence for the presence of a glueball. This has led 
to several phenomenological attempts [16, 21] to describe these three observed 
states in terms of the lattice input. 



8 



As we emphasised above, in fuU QCD on a lattice one just obtains values for 
the ao and /o masses. In the simplified case oí Nf = 2 flavours of degeneráte 
quark, one docs indccd find [9] two /o statcs, and they can be interpreted as 
mixtures of the qq and glueball states with the qq statě having the properties 
found for the ao- 

One useful lattice input would be a dctermination of the oq mass as the quark 
mass is varied in fuU QCD, especiahy because of the problems of determining the 
ao mass in the quenched approximation. At present, the fuU QCD studies [7, 9] 
are hmited to relatively coarse lattice spacing, so the continuum limit is not 
close. Furthermore, as the quark mass is reduced the ao can decay (to nr]) and 
this will influence the lattice analysis [24]. 

3 Hybrid Mesons 

A hybrid meson is a meson in which the gluonic degrees of freedom are excited 
non-trivially. The most dircct sign of this would bc a spin-cxotic meson, since 
that could not be created from a qq statě with unexcited glue. A spin-exotic 
meson could, however, be a qqqq or meson- meson statě and that possibility will 
be discussed. Wc first discuss hybrid mesons with static hcavy quarks where 
the description can be thought of as an excited colour string. The situation 
concerning light quark hybrid mesons is then summarised 

3.1 Heavy quark hybrid mesons 

Considcr QQ states with static quarks in which the gluonic contribution may be 
excited. We classify the gluonic fields according to the symmetries of the systém. 
This discussion is very similar to the description of electron wave functions in 
diatomic molecules. The symmetries are (i) rotation around the separation axis 
z with representations labelled by Jz (ii) CP with representations labelled by 
g{+) and u{—) and (iii) CTZ. Here C interchanges Q and Q, P is parity and TZ 
is a rotation of 180° about the mid-point around the y axis. The CTZ operation 
is only relevant to classify states with J^, = 0. The convention is to label states 
of Jz = 0, 1, 2 by S, n, A respectively. The ground statě will have Jz = O 
and CP = +. 

The exploration of the energy levels of other representations has a long 
history in lattice studies [25, 27]. The first excited statě is found to be the IIu. 
This can be visualised as the symmetry of a string bowed out in the x direction 
minus the same deflection in the —x direction (plus another component of the 
two-dimensional representation with the transverse direction x replaced by y), 
corresponding to flux states from a lattice operátor which is the difíerence of 
U-shaped paths from quark to antiquark of the form Fl — U. 

The picture of the gluon flux between the static quarks suggcsts that the 
excited states of this string may approximate the excited potentials found from 
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Figuře 3: The potential energy between static quarks at separation i? (in units 
of ro « 0.5 fm) [27]. The symmetric gluonic field configuration is shown by 
thc lowcr points whilc thc n^, cxcitcd gluonic configuration is shown above. 
The energy levels in these potentials for 6 quarks are shown using the adiabatic 
approxiniation. 



thc latticc. In the simplest string thcory, thc first cxcitcd Icvcl has n„ symmctry 
and is at energy tt/í? above thc ground statě. This is indeed approximately 
valid and a closer approximation is to use a relativistic version [26] (námely 
E^(R) = (cr2i?2 + 27rcr(m - l/12))i/2 for the m-th level), see also ref. [28] for a 
recent comparison of this expression. 

Recent lattice studies [28] have used an asymmetric space/time spacing 
which cnablcs cxcitcd statcs to bc dctcrmincd comprchcnsivcly. These rcsults 
confirm the finding that the n„ excitation is the lowest lying and hence of most 
relevance to spectroscopy. 

From the potential corresponding to these cxcitcd gluonic statcs, onc can 
determine the spectrum of hybrid quarkonia using the Schrodinger equation in 
the Born-Oppenheimer approximation. This approximation will be good if the 
heavy quarks move very little in thc timc it takcs for thc potential between 
them to become established. More quantitatively, we require that the potential 
energy of gluonic excitation is much larger than the typical energy of orbital or 
radiál excitation. This is indeed the case [25], especially for 6 quarks. Another 
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nice feature of this approach is that thc sclf energy of thc static sourccs cancels 
in thc cncrgy difference between this hybrid statc and thc QQ statcs. Tlius the 
latticc approach gives directly the excitation energy of each gluonic excitation. 

Thc n„ symmetry statě corresponds to excitations of thc gluonic ficld in 
quarkonium caUed magnetic (with L^'~" — 1^ ) and pseudo-clcctric (with 1 ^) 
in contrast to the usual P-wave orbital excitation which has L^'~' = 1 . Thus 
we expect different quantum number assignments from those of the ghionic 
ground statě. Indecd combining with the heavy quark spins, wc get a degeneráte 
set of 8 States: 

Note that of these, J^^ = 1"+, 0+" and 2+- are spin-exotic and hence will 
not mix with QQ states. They thus form a very attractive goal for experimental 
searches for hybrid mesons. 

The eightfold degeneracy of the static approach will be broken by various 
corrections. As an example, one of the eight degeneráte hybrid states is a pscu- 
doscalar with the heavy quarks in a spin triplet. This has the same overall 
quantum numbers as the S-wave QQ statě (r/fc) which, however, has the heavy 
quarks in a spin singlet. So any mixing between these states must be medi- 
ated by spin dependent interactions. These spin depcndcnt interactions will be 
smaller for heavier quarks. It is of interest to establish the strength of these 
effects for b and c quarks. Another topič of interest is the splitting between the 
spin exotic hybrids which will come from the different energies of the magnetic 
and pseudo-clcctric gluonic excitations. 

One way to go beyond the static approach is to use the NRQCD approx- 
imation which then enablcs the spin dependent effects to be explored. One 
study [28] finds that the L^*^ = 1"*"" and 1 excitations have no statistically 
significant splitting although the 1"' excitation does lie a little lighter. This 
would imply, after adding in heavy quark spin, that the = 1 hybrid 

was the lightest spin exotic. Also a relatively large spin splitting was found [29] 
among the triplet states considering, however, only magnetic gluonic excita- 
tions. Another study [36] explores the mixing of non spin-exotic hybrids with 
regular quarkonium states via a spin-flip interaction using lattice NRQCD. 

Confirmation of the ordering of the spin exotic states also comes from lattice 
studies with propagating quarks [30, 31, 34] which are able to measure masses 
for all 8 states. We discuss that evidence in more detail below. 

Because of the similarity of the lightest hybrid wavefunction with that of the 
2S State (which has L = 1), it is convenient to quote mass differences between 
these states. Within the quenched approximation, the lattice evidence for bb 
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Figuře 4: The potential energy for quenched and 2 flavours of sea quark for the 
ground statě and first excited gluonic statě [6]. 

quarks points to a lightcst hybrid spin exotic with J^'-' = 1 ^ at an energy 
given by {niH ~ iti2s)'''q —1.8 (static potential [27]); 1.9 (static potential [28], 
NRQCD [29]); 2.0 (NRQCD [28]). These results can be summarised as {tuh - 
Tn2s)fo = 1-9 ± 0.1. Using the experimental mass of the T(2S'), this implies 
that the lightest spin exotic hybrid is at mn — 10.73(7) GeV including a 10% 
scale error. Above this energy there will be many more hybrid states, many of 
which will be spin exotic. 

The results from a study with N f = 2 flavours of sea-quarks show very little 
change in the static potential (sec fig. 4) and relatively little change in NRQCD 
determinations [29] of mass ratios such as {mu — ^23) / {mip — mis). Expressed 
in terms of ro (using ro = 1.18/^/ct) this gives {mn — m2s)f'o — 2.4(2), however. 
This is significantly larger than the quenched result and, using the IP— 15* mass 
difference to set the scale, yields a prediction [29] for the lightest hybrid mass 
of 11.02(18) GeV. 
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3.2 Hybrid meson decays 

Within this static quark framework, one can explore the decay mechanisms. One 
speciál feature is that the symmetries of the quark and colour fields about the 
static quarks must bc prcscrvcd cxactly in dccay. This has the conscqucncc that 
the decay from a n„ hybrid State to the open-ř» mesons (-BŠ, B*B, BB*, B*B*) 
will be forbidden [41] if the hght quarks in the B and B* mesons are in an S- 
wavc rclative to the hcavy quark (sincc the hnal statě wiU havc the hght quarks 
in either a triplet with the wrong CP or a singlet with the wrong where z is 
the interquark axis). The decay to B**-mesons with light quarks in a P-wave is 
allowed by symmctry but not energctically. 

The only allowed decays are when the hybrid statě de-exites to a non-hybrid 
statě with the emission of a light quark-antiquark pair. Since the n„ hybrid 
statě has the heavy quark-antiquark in a triplet P-wavc statc, the resulting 
non-hybrid statě must also be in a triplet P-wave since the heavy quarks do 
not change their statě in the limit of very heavy quarks. Thus the decay for h 
quarks will be to Xb + M where M is a light quark-antiquark meson in a flavour 
singlet. This proceeds by a disconnected light quark diagram and it would 
be expected [40] that the scalar or pseudoscalar meson channels are the most 
important (ie they have the largest relative OZI-rule violating contributions). 
Lattice estimates [41] of these transitions have been made and the dominant 
mode (with a width of around 100 MeV) is found to be with M as a scalar 
meson, námely ií — > + /o- 

These estimates are in the static quark limit, in which the spin-exotic and 
non spin-exotic hybrid mesons are degeneráte. For the latter, however, the 
interpretation of any observed states is less clear cut, since they could be con- 
ventional quark antiquark states. Moreover, the non spin-exotic hybrid mesons 
can mix directly (ie without emission of any meson M) with conventional quark 
antiquark states oncc one takcs into account corrcctions (of order 1/Mq) to the 
static approximation applicable for heavy quarks with physical masses. 

It is encouraging that the decay width comes out as relatively small, so 
that the spin-exotic hybrid states should show up experimentally as sufficiently 
narrow resonances to be detectable. This decay analysis does not take into 
account heavy quark motion or spin-flip and these effects will be significantly 
more important for charm quarks than for 6-quarks. 

3.3 Light quark hybrid mesons 

I now focus on lattice results for hybrid mesons made from light quarks using 
fuUy relativistic propagating quarks. There will bc no mixing with qq mesons 
for spin-exotic hybrid mesons and these are of speciál interest. The first study of 
this area was by the UKQCD Collaboration [30] who used operators motivated 
by the heavy quark studies referred to above to study all 8 J^^ values coming 
from L^^ = 1+" and 1~+ excitations. The resulting mass spectrum gives 
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the J^'~^ = 1~+ statě as the lightest spin-exotic statě. Taking account of the 
systematic scale errors in the lattice determination, a mass of 2000(200) MeV 
is quoted for this hybrid meson with sš Hght quarks. Although not dircctly 
measured, the corresponding light quark hybrid meson would be expected to be 
around 120 MeV hghter. 

A second lattice group has also cvahiatcd hybrid meson spcctra with propa- 
gating quarks from quenched latticcs. Thcy obtain [31] masses of the 1 statě 
with statistical and various different systematic errors of 1970(90) (300) MeV, 
2170(80) (100) (100) MeV and 4390(80) (200) McV for nň, sš and cč quarks re- 
spectively. For the 0^ spin-exotic statě they have a noisier signál but evidence 
that it is heavier. They also explore mixing matrix elements between spin-exotic 
hybrid states and 4 quark operators. 

The first analysis [34] to determine the hybrid meson spectrum using full 
QCD used Wilson quarks. The sea quarks used had several different masses and 
an extrapolation was made to the limit of physical sea quark masses, yielding a 
mass of 1.9(2) GeV for the lightest spin-exotic hybrid meson, which again was 
found to be the 1 '". In principle this calculation should take account of sea 
quark effects such as the mixing between such a hybrid meson and qqqq states 
such as r]TT, although it is possible that the sea quark masses used are not light 
enough to explore these features. 

A rcccnt dynamical quark study from 2-1-1 flavours of improved staggered 
quarks has also produced results [32]. They also compare their results with 
quenched calculations and find no significant diíference, except that the am- 
bigTiity in fixing the lattice energy scale is better controUed in the dynamical 
simulation since different reference observablcs are closer to experiment. Their 
summary result for the 1 ^ hybrid with strange quarks is 2100 ± 120 MeV, in 
agreement with carlier results. They note that the energies of two-meson states 
(such as TT -|- 6i OT K + ^^(1"*") ) with the hybrid meson quantum numbers are 
close to the energies they obtain. This suggcsts that these two-particle states, 
which are allowcd to mix in a dynamical quark trcatmcnt, may be influencing 
the masses dctermined. A study of hybrid meson transitions to two particle 
states is necded to illuminate this area, using techniques such as those used for 
hcavy quark hybrid decay [41] and decays of light quark vector mesons [42]. 

The lattice calculations [30, 31, 34, 32, 35] of the light hybrid spectrum are 
in good agreement with each other. They imply that the natural energy range 
for spin-exotic hybrid mesons is around 1.9 GeV. The J^^ = 1"+ statě is found 
to be lightest. It is not easy to reconcile these lattice results with experimental 
indications [37] for resonances at 1.4 GeV and 1.6 GeV, especially the lower 
mass value. Mixing with qqqq states such as í^tt is not included for realistic 
quark masses in the lattice calculations. Such effects of pion loops (both real and 
virtual) have been estimated in chiral perturbation theory based models [38] and 
they could potentially reconcile some of the discrepancy between lattice mass 
estimates (with light quarks which are too heavy) and those from experiment. 
This can be interpreted, dependent on one's viewpoint, as either that the lattice 
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calculations are incomplete or as an indication that the experimental states may 
have an important meson-meson component in them. 

The light quark technique of using relativistic propagating quarks can also be 
extended to charm quarks, as was note above [31]. Another group has explored 
the charm quark hybrid states also using a fully relativistic action, albeit with an 
anisotropic latticc formulation [33]. Their quenched study is in agrcement with 
the isotropic lattice result quoted above, finding a mass value of 4.428(41) GeV 
in the continuum limit for the 1 hybrid where the scale is set by the ^Pi — lS 
mass splitting (458.2 McV expcrimcntally) in charmonium. Thcir result is also 
consistent with that from NRQCD methods [29] applied to this case. These 
results all have the usual caveat that in quenched evaluations the overall mass 
scale of the energy difference from the 15 statě at 3.067 GeV is uncertain to 
10% or SO (for example the {2S-1S)/CPi - IS) is found to be 15% higher than 
experiment) which is a major source of systematic error (approximately ±140 
MeV). They also produce estimatcs for other charmonium spin-exotic states: 
0+- at 4.70(17) GeV and 2+" at 4.89(9) GeV. The 0~ statě is not resolved. 

Thus masses near 4.4GeV are found for the charmonium 1 statě using 
relativistic quarks. The non-relativistic approach using NRQCD is expected to 
have big systematic errors for quarks as light as charm, but results [29] do agree 
with this value. The heavy quark eífective theory approach has a leading term 
which corresponds to a static heavy quark, resulting in an estimate [27] of the 
spin-exotic charm statě mass of 4.0 GeV. Here again the systematic error is 
potentially large for charm quarks. 

4 Hadronic molecules 

By exotic statě we mean any statě which is not dominantly a qq or qqq statě. 
For example, a statc madc from hadrons bound in a molcculc would be exotic. 
Examples of hadronic molecules have been known for a long time: the deuteron 
is a proton-neutron molecule for example. It is very weakly bound (2 MeV) and 
is quitc extended. It is more efBciently described in terms to a neutron and a 
proton than as six quarks. 

The residual hadronic interaction, the force between two colour-singlet hadrons, 
is much weaker than the colour force between quarks. Although it is called a 
'strong interaction', it is relatively weak. At large distance, it will be dom- 
inated by the exchange of the lightest hadrons allowed (typically one or two 
pions). For example, the scale of nuclear binding is around 8 McV whcrcas the 
gluonic forces binding three quarks to make a nucleon contribute most of its 
mass of 938 MeV. For this reason lattice methods need to be developed spe- 
cially to tackle this problém. Basically, one is interested in binding energies, so 
it is the energy difference between the two hadrons and the hadronic molecule 
that is of interest. This difference can sometimes be determined better than the 
total energy itself. Even so, the detailed dynamics of such molecular states will 
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depend on the long range forces (typically one or two pion exchange) and this 
will be modified considerably in lattice studies with light quark masses which 
are too hcavy (typically down to 50% of the strange quark mass only). So only 
qualitative input can be obtained from the lattice, but this can still be used to 
validate models. 

Because of the small binding cncrgy and the dominance of pion exchange 
in the binding, it is not feasible to obtain the deuteron binding direct from 
QCD using lattice methods at present. There has been speculation that other 
di-baryon systems might be more strongly bound: the H dibaryon (a AA statě) 
being the best known. If it were strongly bound then it might be stable to 
weak decay which might have astrophysical consequences. The current status 
of lattice studies [43] is that finite box size effects are large but there is no 
convincing lattice evidence that this statě is bound. At the largest volumes 
studied, with L pa 4fm, in quenched simulations the ratio (mír/2mA) — 1 is 
found to be positive and in the range 5 to 15 %. 

Other molecular states involving two hadrons have been conjectured. Several 
meson resonances are known which are dosely connected with nearby thresholds: 
the A(1405) which is just below the KN threshold and the ao(980) and /o(980) 
which are close to the KK threshold. Another statě close to a threshold is the 
A''(1535) which is just above the r]N threshold. Again lattice studies are not 
able to shcd vcry much light directly on these states since the quark masses 
used in the lattice studies are unphysical. However, if they are not produced in 
a lattice study which explores qq and qqq states, this may help to support the 
conclusion that they are primarily molecular in structure. 

One case which is relatively easy to study is the BB systém, idcalised as 
two static quarks and two light quarks. Then a potcntial as a function of the 
separation R between the static quarks can bc dctermincd. Because the static 
quark spin is irrelevant, the states can be classified by the light quark spin and 
isospin. Lattice results [44] (using a light quark mass close to strange) have been 
obtained for the potcntial cncrgy for = 0, 1 and Sq = 0, 1. For vcry heavy 
quarks, a potcntial bclow 2Mb will imply binding of the BB molccules with 
these quantum numbers and L = Q. For the physically rclcvant case of b quarks 
of around 5 GcV, the kinctic cncrgy will not bc ncgligiblc and the binding energy 
of the BB molecular states is less clcar cut. One way to estimate the kinetic 
energy for the BB case with reduced mass circa 2.5 GeV is to use anály tic 
approximations to the potentials foimd. For example the Iq, Sq={0,0) case (see 
fig. 5) shows a deep binding at i? = O which can be approximated as a Coulomb 
potcntial of — 0.1/i? in GeV units. This will give a di-mcson binding energy of 
only 10 MeV. For the other interesting case shown in fig. 6, {Iq, Sq)=(0,í), a 
harmonie oscillator potcntial in the radiál coordinate of form — 0.04[1— (r— 3)^/4] 
in GeV units leads to a kinetic energy which completely cancels the potcntial 
energy minimum, leaving zero binding. This harmonie oscillator approximation 
lies above the estimate of the potcntial, so again we expect weak binding of the 
di-meson systém. 
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Figuře 5: Th(^ binding cncrgy [44] bctwccn two hcavy-light mcsons (with static 
heavy quarks and liglit quarks of mass corrcsponding to strange) at scparation 
R (in units of ro ~ 0.5 fm) with the two light quarks having 1=0 and S=0. 



Bccausc of these vcry smaU valucs for the di-mcson binding cncrgics, one 
needs to retain corrections to the heavy quark approximation to make more 
definite predictions, since these corrections are known to be of magnitude 46 
McV from the B, B* sphtting. It will also be necessary to cxtrapolatc the 
hght quark mass from strange to the hghter u, d values to make more definite 
predictions about the binding of BB molecules. 

Modcls for the binding of two B mcsons involvc, as in the casc of the 
deuteron, pion exchange. The lattice study [44] is able to make a quantita- 
tive comparison of lattice pion exchange with the data described above using 
lattice dctcrminations of the B* Btt coupling [45] and cxcellent agreement is 
obtained at larger R values as shown in fig. 7, as expected. 

5 Conclusions and Outlook 

Quenched lattice QCD is well understood and accurate predictions in the con- 
tinuum limit are increasingly becoming available. The lightest glueball is scalar 

with mass to(0"'"+) — 1611 (30) (160) McV whcre the second crror is an ovcrall 
scale error. The excited glueball spectrum is known too. The quenched ap- 
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Figuře 6: Th(^ binding cncrgy [44] bctwccn two hcavy-light mcsons (with static 
heavy quarks and liglit quarks of mass corrcsponding to strange) at scparation 
R (in units of ro ~ 0.5 fm) with the two light quarks having 1=0 and S=l. 



proximation also givcs Information on quark-antiquark scalar mcsons and thcir 
mixing with glueballs. This determination of the mixing in the quenched ap- 
proximation also sheds Hght on results for the spectrum directly in full QCD 
whcrc the mixing will bc cnablcd. In full QCD, the scalar meson masses are 
determined directly but there is no concept of a glueball as such, much as in 
the experimental case. Additional work is need to reduce the lattice spacing, or 
use improvcd actions, to cxplore the continuum limit for scalar mcsons in full 
QCD. There is also some lattice Information on the hadronic dccay amplitudes 
of glueballs and this is an area where further study may be anticipated. 

For hybrid mcsons, there will bc no mixing with qq for spin-cxotic statcs 
and these are the most useful predictions. The J^'^' = 1 ^ statě is expected 
in the range 10.7 to 11.0 GeV for b quarks, 2.0(2) GeV for s quarks and 1.9(2) 
GcV for u, d quarks. Mixing of spin-cxotic hybrids with qqqq or cquivalcntly 
with meson-meson is allowed and will modify the predictions from the quenched 
approximation. A first lattice study has been made of hybrid meson decays. For 
heavy quarks, the dominant mode is string de-cxcitation to x + /o where /o is a 
flavour singlet scalar meson (or possibly two pions in this statě). The magnitude 
of the decay rate is found to be of order 100 MeV, so this decay mode should 
still leave a detectably narrow resonance to be observed. 
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Figuře 7: The contribution [44] to the binding energy for the spin and isospin 

combinations corresponding to tt cxchangc (octagons) and p cxchangc (fancy 
squares). The sohd hne gives the tt exchange contribution which is normahsed 
by the B*Bn couphng. The p exchange prediction has a free normahsation 
and is shown by the dash-dottcd Hne. The results are plotted versus interquark 
separation R (in units of a « 0.17 fm). 

The topič of possible multi-quark bound states is diíBcult because the scale 

of the cxpcctcd binding cncrgics is a fcw McV and this small valuc is a challcnge 
for lattice studies. As an example, some evidence was presented for a possible 
BgBs molecular statě. 
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